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METHOD OF ANALYSIS FOR COMPRESSIBLE FLOW THROUGH MIXED-FLOW CENTRIFUGAL 

IMPELLERS OF ARBITRARY DESIGN 1 

By Joseph T. Hamrick, Ambrose Ginsburg, and Walter M. Osborn 


SUMMARY 

A method is presented for analysis of the compressible flow 
between the hub and the shroud of mixed-flow impellers of 
arbitrary design. Axial symmetry was assumed, but the forces 
in the meridional ( hub to shroud) plane, which are derived from 
tangential pressure gradients , were taken into account. 

The method was applied to an experimental mixed-flow im- 
peller. The analysis of the flow in the meridional plane of 
the impeller showed that the rotational forces, the blade curva- 
ture, and the hub-shroud profile can introduce severe velocity 
gradients along the hub and the shroud surfaces. Choked 
flow at the impeller inlet as determined by the analysis was 
verified by experimental results. 

INTRODUCTION 

Because of the complexity of the numerical solution for a 
three-dimensional analysis of flow through a centrifugal com- 
pressor, most treatments of the problem have been concerned 
with a two-dimensional analysis of the flow between the 
pressure side of one impeller blade and the suction side of 
the preceding blade. Relaxation methods in which it was 
assumed that no velocity variation existed between the hub 
and the shroud have been used (reference 1) to give an 
analysis of the flow between blades. In the conventional 
mixed-flow centrifugal impeller, however, the nonsymmetri- 
caf hub-shroud profile, the nonuniform energy addition 
across the passage, and the blade-angle variation due to close 
adherence to radial blade elements may introduce velocity 
gradients between the hub and the shroud comparable with 
those between the blades. For mixed-flow compressors, 
therefore, an analysis of the flow in the meridional plane is 
required. (The meridional plane is one that passes through 
and is parallel to the axis of rotation.) 

In order to investigate the isentropie compressible flow in 
the meridional plane, a flow equation based on equilibrium, 
continuity, and general energy considerations was derived 
at the XACA Lewis laboratory during 1949-50 and is 
presented herein. Axial symmetry was assumed, but the 
forces that result from tangential pressure gradients were 
taken into account. The numerical solution involves the 
adaptation of stream-filament theory and the numerical 
integration of all equations. 

The applicability of the method for predicting the flow 
witliin the impoller was investigated by analyzing the flow 


through a mixed-flow centrifugal impeller for which experi- 
mental data are available (reference 2). Flow streamlines, 
pressure distributions, and velocity distributions relative to 
the impeller for a given set of operating conditions are 
presented. Theoretical shroud pressures and theoretical 
maximum weight flows over a range of speeds are compared 
with those experimentally obtained. 

ANALYSIS 

An analysis is made in the meridional plane of the mixed- 
flow impeller. Axial symmetry is assumed, but forces in the 
meridional plane that are derived from tangential pressure 
gradients are taken into account. By the use of the methods 
presented in this analysis, the isentropie compressible flow 
in the meridional plane may be obtained and shown in 
terms of the flow streamlines and the velocity and pressure 
distributions. 






(a) Impeller profile in meridional plane. 

(b) Intersection of developed cone surface. 

Figure I —View of impeller showing orientation of blade angle and forces acting on rotating 
fluid particle in meridional plane. 


1 Supersedes NACA TN 2165, “Method of Analysis for Compressible Flow Through Mixed- Flow Centrifugal Impellers of Arbitrary Design’' by Joseph T. Hamrick, Ambrose Giosburg, 
and Walter M. Osborn, 1950, 
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Geometrical quantities. — A complete geometrical descrip- 
tion of the blade is necessary for a numerical solution. 
Geometrical quantities are defined in appendix A and shown 
in figures 1 and 2. The blade angle /3 is oriented in the 
following manner: In figure 1(a), a line OB, which intersects 
the axis of rotation at 0, is drawn tangent to the point P on 
a streamline in the meridional plane. If the tangent line 
OB is rotated about the center line of impeller rotation O-O' , 
a conic surface of revolution that intersects the impeller 
blades is formed. This intersection is shown on a developed 
surface of the cone in figure 1(b). At the point of tangeney 
to the streamline (point P), the blade mean line forms the 
angle — /3 with the tangent line- 

limitations and assumptions.- — The fluid is assumed to be 
nonviscous but compressible. All compression is by the 
adiabatic, isentropic, steady-flow process. The velocity is 
assumed uniform from blade to blade in the direction of 
rotation (axial symmetry) and the average angle of flow 
relative to the impeller is assumed the same as the blade 
angle jS. This assumption is considered valid where channel 
flow prevails, but is questionable at the inlet and the outlet. 



W 





(a) Meridional plane components. 

(b) Components of force normal to blade. 

(c) Orientation of <p with respect to £. (Dashed lines lie in meridional plane.) 

Figure 2 —Quantities used in development of force equations. 


Method of analysis. — The method is based on one proposed 
by Fliigel for analysis of flow in the blade-to-blade direction 
and discussed by Stodola in reference 3. The forces acting 
on a particle along the streamline path are considered. The 
component of the blade surface force normal to the meri- 
dional streamline, the centrifugal forces resulting from the 


curvilinear path of the particle, and the force due to rotation 
about the axis of the impeller are equated to the change in 
pressure across the particle in the hub-to-shroud direction. 
The resulting force equation is combined with one based on 
the general energy equation to give an equation involving 
the relative velocity q and a number of variables dependent 
upon the geometry of the impeller. The variables that lend 
themselves to graphical treatment are lumped into two para- 
meters giving a first-order differential equation linear in q. 
This equation is used in conjunction with the continuity 
equation to obtain the velocity distribution. Streamlines 
between the hub and the shroud similar to those shown in 
figure 1 (a) are drawn in by estimate. When the differential 
equation is solved for q and a velocity is assumed at the hub, 
the velocity distribution from hub to shroud may be obtained 
by numerical integration. By means of the continuity equa- 
tion, the weight flow corresponding to the assumed hub 
velocity is obtained. Several solutions may be necessary for 
obtaining the desired weight flow. The streamlines may be 
rearranged to match the new values of velocity and density 
and the complete analysis repeated. 

Force equation, — For the forces acting on the fluid 
particle in the n direction on the meridional plane (fig. 1 (a)), 
the following equation is obtained (appendix B): 


1 dp v 2 . (wr+u) 2 dr 
p dn~ r c r dn 


From figure 1 it can be seen that 


and 


v = q cos 
u = q sin |3 


( 1 ) 


In terms of q, equation (1) becomes 


dp 
> dn 


,/ cos 2 (3 , sin 2 (3 dr \ . „ dr . dr „ . . 

'Tt+ts <2) 


Energy addition.- — From general energy considerations, 
the isentropic adiabatic energy addition is 

^ [wrAcorj+tt!)— «X,] = Ai ~ [(ediy+u,) 2 -!- !)! 2 ]— H { 

where the prerotation term X=ri(wr i -j-u < ) 


i why 2 , u, 2 , ty 2 TT wX x 

h '-^f + - e + Tr a '-— 


( 3 ) 


From figure 1, 


q 2 —u 2 -\-v 2 


and with the assumption that II t is constant for all stream- 
lines, 

l(m) 2 -tf]+~r = h 0 -j- g [(wr a ) 2 — g 0 2 l + '~ 


hx — hn — ^ [w 2 (r t 2 r o 2 ) (q 2 q o 2 ) 2w(X t Xo)] 
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For an isentropic process 

( n i 
— - dp 
opg 

Therefore 

C 1 - dv=iz [«*( r ffo 5 ) — 2 w(X L — Xq)] (4) 

Jo Pi *g 


When differentiated with respect to n, equation (4) becomes 


1 dp „ dr da d\ 
p dn <l3 ~ r ~3n ® dn 03 dn 


( 5 ) 


General velocity equation. — Equations (2) and (5) are 
combined to give the following equation: 



/ cos 1 
\ r c 


sin 1 0 dr\ 
r dn) 


dr F 

-2w sinjS ■ - 

d n q 


u d\ 
q dn 


( 6 ) 


For straight blades lying on the meridional plane with no 

F o) d\ 

prerotation of the fluid at the inlet, the terms — and - -j— 
r ’ q q dn 

drop out; a first-order differential equation linear in q thus 

remains. In order to arrive at a solution for curved-blade 

impellers and for prerotation at the inlet, it may be desirable 

to solve for q and to neglect these terms in the first solution. 

The values of q thus found may be used in the force and pre- 

rotation terms for the second solution and may be treated in 

the same manner as those terms dependent upon the geometry 

of the impeller. 


Velocity equation for impeller of arbitrary blade shape 
having radial blade elements with no prerotation. — The 
blade-force component F K normal to the meridional stream- 
lines is evaluated in appendix C as 

^=— sin 0 tan a ^2<o sin sin* 0 sin a tan a 


The velocity equation for a radial-bladed impeller without 
prerotation becomes 


dq /cos 1 0 sin 1 0 dr sin 1 0 . 

1n*\rz ~ Sl 


sin a tan 




2u sin 0 sin 0 tan a (2 a sin (7) 


Variables dependent upon the impeller geometry or upon 
values obtained from a previous solution are represented by 
parameters a and b where 


cos 1 0 sin 1 0 dr sin 1 0 . 
a—— — -— , sin a tan a 


r c 


r dn 


and 


6= 2w sin 0 ^--|rsin 0 tan a ^2« sin 


Equation (7) becomes 


^S=aq—b 


dn 


( 8 ) 


When equation (8) is solved for q, 

f*oia/ pa —f'adn \ 

q=e^° {g 0 — I be “™ dnj W 

In solving for q, the term du/dm in parameter b is neglected 
for the initial solution in order to obtain estimated values of 
duldm for the succeeding solution. 

Blade surface velocity distribution. — If a linear pressure 
gradient from blade to blade is assumed, the blade surface 
velocity distribution may be approximated for a finite number 
of blades by the method given in appendix D. 

Continuity equation. — The equation for continuity of 
flow through the hub-to-shroud passage may be written 

to=2ffJ* rfpgq cos 0 dn (10) 

where / is a factor that takes into consideration the area 
taken up by the blade cross section. The value of the 
density p is determined from the energy equation (equation 
(3)) and the relations 

jL—( A Y ” 1 

NUMERICAL EXAMPLE 

IMPELLER EXAMPLE 

A numerical example has been computed for an 18-blade 
mixed-flow centrifugal impeller, the performance of which is 
described in reference 2. The impeller was designed for a 
constant cross-sectional flow area from inlet to outlet. The 
blade curvature is such that the intersection of the blade 
with any circular cylinder having the same center line as the 
impeller produces a parabola on the developed cylinder sur- 
face. Because of adherence to radial blade elements, the 
blade angle 0 increases from hub to shroud. 

The operating conditions for the impeller were as follows: 


Equivalent impeller tip speed, U[- <fe, ft/sec 1331 

Ratio of specific heats, y 1.4 

Prerotation at inlet 0 

Impeller equivalent weight flow, Wy'e/S, Ib/sec 8.743 


PROCEDURES FOB NUMERICAL SOLUTION 

The meridional plane of the impeller together with the 
estimated streamlines and resulting normals is shown in figure 
1(a). Estimated streamlines are drawn in between the hub 
and the shroud; equal mass flow is assumed between the 
annular stream tubes. Lines normal to the stre amlin es are 
drawn in along the flow path. 

Radius of curvature. — The radius of curvature of the 
streamline at any point in the meridional plane r t is obtained 
from 

. [l+CrTP 1 

r,- p 

The slope of each streamline r' at each normal n is plotted 
against axial distance in order to obtain the change in slope r". 
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Numerical integration. — The velocity q is found by 
numerically integrating equation (9) 


? 



( 9 ) 


The velocity q 0 is the estimated velocity at the hub. 

When a is plotted against n (fig. 3(a)) where n is the dis- 
tance along the normal from hub to shroud, the integral 

a dn can be evaluated by the trapezoidal rule. A plot of 



-Pain 

be against, n can then be obtained (fig. 3(b)) and 

evaluated in a similar manner. As previously stated, the 
term duldm in parameter b is neglected in the first solution. 
The evaluations of these integrals allow the determination of 
the velocity q at any point n. 

The weight flow at any normal may be obtained from 
equation (10) 


w—2 IT 



cos /3 dn 


( 10 ) 


When a plot of rfpgq cos against n is used (fig. 3(c)), the 

integral I rfpgq cos 0 dn can be evaluated by the trapezoidal 

rule. If the weight flow does not equal the designated 
weight flow for the passage, a new value of the velocity 
is assumed and the solution is repeated. 

In order to satisfy continuity for each stream tube, new 
streamlines are drawn based on the velocities of the preceding 
solution. Velocity q (fig. 4(a)) and density pg (fig. 4(b)) 
are plotted against n. The radii at which the new stream- 
lines cross the normals are obtained as follows: From 


and 


An 


r,+i-r a 

cos a 


Aw=2 vrjpgq cos j3 An 


Aw cos ot 
vfpgq cos jS 



where cos a, f, p, q, and cos (3 are average values between 
streamlines based on estimated values of r n+I . After the 
new streamlines are determined, the normals are drawn. 
The graphical solution is then repeated for succeeding sets 
of new streamlines until continuity is satisfied. 

ACCURACY 

The accuracy of a solution depends to a large extent upon 
the accuracy with which the streamline curvature is deter- 
mined. Measurement of other geometrical quantities, of 
course, enters into the accuracy, but their measurement is 
more precise. 

The sudden decrease in flow area (as much as 25 percent 
at the hub) from a point just ahead of the blades at the inlet 
to a point just inside the blades distorts the streamline pic- 
ture. A tedious and detailed solution is necessary for a true 
picture of the distortion in tliis region. Also, the assumption 
that the flow angle is the same as that of the blades at the 
inlet is questionable. A lengthy analysis at the inlet was 
therefore considered unjustifiable. 

The cross-sectional area for the clearance is approximately 
4 percent near the inlet and affects the maximum flow. 
In the solution, the clearance area was treated as though the 
blades extended to the shroud because of continuity con- 
siderations. 



(a) Parameter®. 


- f“ a in 
e Jo 


(b) Exponential be 

(c) Weight -flow factor. 

Fiodei 3.— Plots used far numerical Integration of equations (10) and ( 11 ). 
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<ai Relative velocity. 

<b) Density. 

FV-vbe +.— Plots used for computing new streamlines. 


RESULTS AND DISCUSSION 

A solution of the flow conditions throughout the passage 
from hub to shroud is shown in figure 5. 

Streamlines. — The flow streamlines in the meridional 



Ca) Streamlines for compressible flow. 

Floral S.— Flow analysis for centrifugal Impeller In meridional plane. 


plane are shown in figure 5 (a). Although the solution has 
been made for equal quantities of flow between streamlines, 
little information can be derived from the streamline con- 
figuration, as shown, other than the direction of flow because 
the variation, in distance between streamlines is due to a 
variation in blade thickness and radius as well as in velocity. 
Near the outlet where blade thickness and radius have little 
effect, however, the streamline spacing indicates considerable 
variation in velocity across the passage. 

Velocity ratio. — Lines of constant velocity ratio Q relative 
to the impeller are shown in figure 5 (b). Along the shroud 
(fig. 6), the velocity is near sonic at the inlet (Mach number, 
0.973) with immediate deceleration just inside the impeller 
and continued deceleration to the outlet, a condition that is 
highly conducive to boundary-layer build-up and flow separa- 
tion. A similar condition exists along the straight part of 
the hub (A to B, fig. 5 (b)). Along the curved part of the 
hub (B to C), the velocity variation is small with a slight 
acceleration existing near the impeller outlet. 

The relative-velocity distribution as obtained in the solu- 
tion represents the mean velocity from blade to blade. 
According to the analysis of reference 1, the velocity along 
the suction surface is much larger than that along the pres- 
sure surface; the velocity varies somewhat uniformly across 
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Velocity 
ratio, (, 


.40 .30 .20 C 



t iQUBB 6.— Concluded. Flow analysis 


Ratio of static 
pressure to inlet 
total pressure 



centrifugal impeller in meridional piano. 


the passage. If the blade unloads at the outlet, as it must 
to satisfy the Kutta condition of tangency at the tip, the 
decelerations along the suction surface of this blade are of 
necessity much greater tlian that shown for the average 
velocity (fig. 5(b)). Such adverse velocity gradients could 
easily cause considerable separation along the flow path. 
At the impeller outlet, these flow separations would result 
in poor flow conditions with consequent losses at the diffuser 
inlet. 

Pressure ratio.— -Lines of constant pressure ratio are 
shown on figure 5 (c). The pressure lines approximately 
parallel the normal lines and indicate that static-pressure 
readings on the stationary shroud of this impoller give repre- 
sentative values for the static pressure across the passage 
from hub to shroud. In figure 7, experimental values of 
pressure ratio along the sliroud are compared with the theo- 
retical values obtained in the numerical solution. Pressure 
ratios are shown for the theoretical weight flow and bracket- 
ing experimental weight flows. The experimental curve for 
the higher weight flow represents a point very near the maxi- 
mum flow or flow cut-off line and shows a large pressure loss 
near the inlet. The experimental curve for the lower weight 
flow represents an operating point near maximum efficiency 
and consequently approaches the theoretical conditions. 
No pressure drop occurs at the inlet, which is in agreement 


with the theoretical curve; however, the theoretical pressure 
increases more rapidly along the sliroud because of ineffi- 
ciency due to friction and flow separation, and hecauso of a 
velocity increase resulting from flow-area reduction duo to 
boundary layer and from a slightly liigher theoretical tip 
speed. 

Maximum flow. — The maximum theoretical weight flow 
across any normal can be determined by making a series of 
estimates of velocity at the hub and obtaining the corre- 
sponding weight flows. For the numerical example, tho 
greatest restriction of flow occurred just inside the inlet. 
The maximum flow over a range of tip speeds was deter- 
mined by the consideration of this flow restriction. Inas- 
much as the radius of curvature was infinite for the stream- 
lines near the inlet in the complete numerical solution, it 
was taken as infinite in computing maximum weight flow 
at all speeds. A comparison of the theoretical and experi- 
mental maximum flow over the speed range is shown in 
figure 8 (a). The slightly higher theoretical values are in 
good agreement with the experimental values over a range 
of tip speeds up to 1400 feet per second (deviation of the 
order of 3 percent). At speeds above 1400 feet per second, 
the deviation increased to a maximum value of 7 percent at 
a tip speed of 1600 feet per second. 
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Fiotm* 6.— Comparison at relative velocity and blade angle along shroud. 



Figuki 7— Comparison of theoretical pressor® ratio along shroud with experimental values 
tor varying weight How. 


The high relative velocities within the impeller, which 
caused choking, resulted from the large blade angle (3 near 
the impeller inlet. A large inlet blade angle was required 
for shockless inlet conditions resulting from the considera- 
tion of the impeller mass flow and the rotational speed. 
These high relative velocities, however, can be reduced by 
an appropriate selection of blade shape to give a more rapid 
decrease in /J near the impeller inlet. A comparison of j3 




(») Comparison of maxim nm flow capacity. 

(b) Comparison of blade-angle decrease from Inlet to outlet at shroud. 

Fiocbi 5.— Comparison of change In blade angle and maximum flow capacity of two Impellers 
having same hub-shroud profiles and blade Inlet angles. 

for the example impeller A and a second impeller B (refer- 
ence 3) is shown in figure 8(b); j3 for each impeller is plotted 
against the distance ratio L along the shroud. Impeller B 
has the same hub-to-shroud profile and inlet blade angle (3 
as that of the impeller A ; however, as indicated in figure 8 (b), 
the rate of decrease in blade angle for impeller B is 
much larger near the impeller inlet. The resulting higher 
flow capacity of impeller B (fig. 8 (a)) may be attributed to 
the larger rate of decrease of blade angle and consequently 
lower relative velocities in the impeller passage near the 
inlet. 

SUMMARY OF RESULTS 

A method has been developed for a compressible flow 
analysis between the hub and the shroud of mixed-flow 
impellers of arbitrary design. Axial symmetry was assumed, 
but the forces in the meridional plane that are derived from 
tangential pressure gradients were taken into account. 
Application of the method to a mixed-flow impeller for 
which experimental data are available yielded the following 
results: 

1. A large variation in velocity existed across the passage 
from hub to shroud. 

2. Adverse velocity gradients occurred along both hub 
and shroud surfaces. 
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3. The theoretical pressure ratio increased more rapidly 
than the experimental along the shroud from inlet to outlet. 
Inefficiency of the actual compression process due to flow- 
separation accounts for the greater part of this deviation. 

4. The maximum flow over a range of speeds as theoreti- 
cally determined closely agreed w-itli the experimental. 

5. For the example calculation, values of static pressure 
at any point along the shroud were representative of those 
across the flow path from hub to shroud. 


CONCLUSION 

Analysis of flow in the meridional plane of a mixed-flow 
impeller showed that the rotational forces, the blade curva- 
ture, and the hub-sliroud profile can introduce adverse 
velocity gradients along the hub and shroud surfaces. 

Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, April IS, 1950 
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APPENDIX A 

SYMBOLS 


The follow-ing symbols are used in the calculations and 
the figures: 

. cos 2 5 sin 2 13 dr sin* § . 
parameter equal to — — — — sin a tan a 

• di* 

parameter equal to 2« sin 9 ^-+sin 8 tan a 



speed of sound for stagnation conditions, ft/sec 
force per unit mass 

factor accounting for area taken up by blade cross 
section 

acceleration due to gravity, ft/sec* 
absolute total enthalpy, Btu/lb 
static enthalpy, Btu/lb 
mechanical equivalent of heat 

ratio of distance m at any point to m at impeller exit 
distance along streamlino from inlet, ft 
number of blades 

distance along normal to streamlino in meridional plane, 
ft 

total pressure, Ib/sq ft 
static pressure, lb/sq ft 
velocity ratio, g/c< 

velocity relative to impeller (fig. 1(b)), ft/sec 
radius of particle from axis of impeller rotation, ft 
radius of curvature of streamlines in meridional plane 
(positive in fig. 2(a)), ft 

slope of streamline in meridional plane with respect to 
axis of rotation - - 

change in slope r' with respect to axis of rotation 
time, sec 


U actual impeller tip speed, ft/sec 

u tangential velocity relativo to impeller (positive in 
direction of rotation, fig. 2(c)), ft/sec 
v through-flow- component of velocity (fig. 2(c)), ft/sec 
W total compressor flow- rate, lb/sec 
w flow rate, Ib/sec 
z distance in axial direction, ft 

a angle betw-een tangent to streamline in meridional 
plane and axis of rotation, deg 
8 blade angle (fig. 1, negative for backward-sw-ept blades), 
deg 

y ratio of specific heats 

5 ratio of total pressure at inlet to standard sea-level 
pressure 

8 ratio of inlet stagnation temperature to standard sea- 
level temperature 

X prerotation term equal to r i (u>r { + u t ) 
p mass density, lb-sec 2 /ft 4 

p c total mass density, lb-scc s /ft 4 

<p blade anglo in axial-tangential direction (negative for 
backward-swept blades, fig. 2(c)), deg 
« angular velocity of impeller, radians/scc 

Subscripts: 
a average 

d pressure or driving face of blade 

i inlet conditions 

n indicates position along normal from hub to shroud or 
normal component 
r radial component 

8 suction or trailing face of blade 

z axial component 

9 tangential component 

0,1 points at succeeding streamlines along normal 


APPENDIX B 


DEVELOPMENT OF FORCE EQUATION IN HUB-TO-SHROUD PLANE FOR FLOW WHEN AXIAL SYMMETRY IS ASSUMED 

Lorenz’s equations of motion relative to the rotating 
impeller, in which axial symmetry is assumed (reference 3), 
are 


F _A dp I (cor+tt ) 2 

r p br^dt* r 

* pbz^dt 2 
v du . u dr dr 

Fe= lt + TTt +2a di 


where the blade force F is the force per unit mass. 
From figure 2(a), 


F n =F, sin a— F t cos a 

where F a may be evaluated from Ft (appendix C) as 

d i z . 1 bp d i r 

dt 2 

(tdr-f-u) 8 


1 bp . 
p b< 


F„=— ~ sin a+' :: j 7 r sin a cos a- 


p br 


dt 1 


COS a+ 


COS a 


(Bl) 
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dr 

v sin a 



Substituting equations (Bla) to (Blc) in 

dt 




dz_ 

dt 

v cos a 



„ 1 dp da , (ar+u)' 

pdn ° dt 1 r 

d l r da . . 

T-,=r cos a -77+sm 
at~ dt 

d v 
a lt 

(Bla) 

From figure 2(a), 

, v dt 

ua= 

r c 

d 2 z 

-ts = — v sm 

at- 

da , dr 

a lt~ COSOL li 

(Bib) 

Equation (B2) therefore reduces to 

dp dz ^_bp dr 

bp . 

1 dp 

COS a 
dr 

(Blc) 

1 dp v 2 , (cor-f-u) 2 dr 

bz dn 1 dr du 

— — sm 
bz 

p dn r c 1 r dn 


COS a 


-F. 


(B2) 


(B3) 


APPENDIX C 


PROCEDURE FOR EVALUATING BLADE FORCE FOR IMPELLER OF ARBITRARY BLADE SHAPE HAVING RADIAL BLADE 

ELEMENTS 


For an impeller with, radial blade elements, a true picture 
of the blade curvature lies on a circular cylinder having the 
same center line as the impeller. The blade force F normal 
to the blade (fig. 2(b)) has the component Ft, which may be 
determined from time rate of change of angular momentum 
in the direction of rotation. The component Ft as given in 
appendix B makes the angle <p with the force F. 

From figure 2(c), 


and 


giving 


F n ——Ft tan <p sin a 

u tan B 

tan <c= = — — 

i: cos a cos a 

F n =—Fe tan j3 tan a 


From appendix B, 



u dr~ , 

7 dt + 


2 


dr 
" dt 


Substituting dm=o dt and drjdm= sin a gives 

Ft=v ( sin a+2co sin a) 

\d m r j 

F n = — v tan jS tan a sin a72w sin 

F n = — o tan jS tan a ^2 to sin si n a tan “ tan £ 

From figure 1, 

r=q cos j9 


u=q sin /5 

F n = — q sin (3 tan a ^2 to sin “7^^— y sin 2 jS sin a tan 


APPENDIX D 

METHOD OF EVALUATING BLADE SURFACE VELOCITIES 

From appendix C, 

„ / du . u . , „ - \ 

F»=v\ s — — sin c(t2(i) sm a 
V dm r ) 


where Ft is the blade force per unit mass in the tangential 
direction. Substituting v—q cos (3. u—q sin (3, and dm= 


- — gives 

cos a 

„ „ t du . q sin |S sin a . _ . \ 

Ft=q cos cos “7* i“2cj sin aj (Dl) 

The quantity Ftp is the blade force per unit volume, which 
when multiplied by the circumferential distance between the 
blades gives the blade force per unit area or the pressure 
difference between the blades; thus, 

F,p(^)=Ap (D2) 

Substituting equation (Bl) in equation (B2) gives 


„ /2? -rf\/du . q sin (3 sin a , „ . \ 

A p = pq cos j8( y AoF C ° S “ ‘ r 72 co sm a) 

(D3) 


Assuming a linear variation of pressure from blade to blade 
yields 

Ap 

p,=Vz — f 


Vi = Va 7nf 


Divi ding by the total pressure P t upstream of the impeller 
inlet gives 

Pt Pa Ap n'li'i 

P~Pt 2 P £ 1 ' 

and 

Pd Pa , Ap 

Pi Pi ' -Pi 

Equation (11) is 

.. . fi i y if 


(Do) 


[7“ r Y 

/if 

> 

Lw 

U.i _ 

gJH t ] 


1 

7-1 
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For zero prerotation at the inlet (X=0) and with 



the expression for velocity becomes 
Z 2 =(ur) 


uecomes 


(DO) 


Substituting equations (B4) and (B5) in equation (B6) gives 
the final equations for determining the velocity on the blade 
surfaces, 



and 
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